Equivalent Hamiltonians 



O 

(N 

=3 



< 



43 



> 



W. N. Polyzou* 

Department of Physics and Astronomy, The University of Iowa, Iowa City, IA 52242 

(Dated: September 1, 2010) 

I give a characterization of the conditions for two Hamiltonians to be equivalent, discuss the 
construction of the operators that relate equivalent Hamiltonians, and introduce variational methods 
that can select Hamiltonians with desirable features from the space of equivalent Hamiltonians. 

PACS numbers: 21.45.Ff,13.75.Cs,03.65.Nk 



I. INTRODUCTION 



In this paper I study the freedom available to redefine interactions without changing the scattering and bound-state 
observables of a quantum mechanical system. I consider how this freedom can be used to formulate interactions that 
have advantages in various situations. The Hamiltonians discussed in this paper are equivalent in the mathematical 
sense at all energy scales; this a stronger requirement than the more flexible notion of equivalence used in effective field 
theory that only requires equivalence up to some order in the expansion parameter. Even with this more restrictive 
fH ' notion of equivalence there is a very large class of equivalent Hamiltonians. 

fS ,. There are a number of formal methods that start from a set of high-quality two and three-body interactions and 
construct an equivalent set of interactions that fit the same bound-state and scattering data. These include renormal- 
ization group methods and methods based on specific unitary transformations that block diagonalize Hamiltonians 
[H H • The compelling feature of all of these methods is that the off-diagonal matrix elements that cou- 
ple the high- and low-energy parts of the problem are suppressed in the transformed interactions. This leads to a 
low-energy effective theories that are approximately decoupled from the high-energy part of the problem. This has 
computational advantages in many-body calculations. The price paid is that the transformed Hamiltonian has new 
many-body forces involving any number of particles. This is similar to what is observed using field redefinitions in 
(yT) ' effective field theories, although the transformed theories discussed in this paper are in principle equivalent to the 
£Nl ! original theories for all energies. 

lf~) • In this paper I introduce a method that can be used to provide independent control of the two, three, and many- 
body interactions. Much of the work contained in this paper has been discussed in Q. The approach is to start 
with the general class of equivalent interactions. This is then restricted to a subset that can be treated variationally. 
Positive functionals are introduced that have minimum values for equivalent potentials with selected properties. For 
. example it is possible to design functionals that select models where the dynamics for energies above some given scale 
approximately decouples from the dynamics for energies below some scale, models that have weak three and four-body 
interactions or models that emphasize an approximate symmetry. 



II. MULTICHANNEL SCATTERING THEORY 

In this section I give a brief summary of multichannel scattering theory that is relevant for this work. 
The Hilbcrt space, H\, for a single particle of mass m and spin j is the space of square intcgrablc functions of 
particle's linear momentum and magnetic quantum number 

(p, = V>(p, n) =[d P J2 IV^P, M)| 2 < oo. (2.1) 

The iV-particlc Hilbcrt space is the AT-fold tensor product of single-particle Hilbert spaces 

H := ®? =l Hi. (2.2) 
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The total linear momentum and total Galilean mass of the iV-particle system are the multiplication operators 

N N 



P :=J2p* M = Y j m i . (2.3) 

i=l i=l 

The iV-body Hamiltonian H has a iV-body bound state if the center-of-mass Hamiltonian, 

has a discrete eigenvalue, — e. If the Hamiltonian is rotationally and translationally invariant it is possible to find 
simultaneous eigenstates of h, p, the total iV-body spin and the projection of the total iV-body spin on the 3-axis. I 
denote these eigenstates by |(ej, ji), p, fi), where the index i labels different bound states when h has more than one 
bound state. 

To define scattering channels let a denote a partition of the N- particles into n a disjoint non-empty clusters of n 0i 
particles. There is a scattering channel a associated with the partition a if there is a n Qi -body bound state in each of 
the n a clusters of the partition a. Channel states asymptotically look like a collection of n a mutually non-interacting 
bound clusters. 

The direct product of the n a bound states in the channel a, 

= |(ei,ii),Pi,Mi) x •■■ x |(e„ a , j„J,p„ a ,/in a ) (2.5) 

defines the mapping $ ai , called the channel injection operator, from the channel Hilbcrt space, H ai , which is tensor 
product of n a single particle Hilbert spaces, 

n ai =Hi®---®H na (2.6) 

to the iV-particle Hilbcrt space by 

, /il) X • • ■ X | {h na , jn a ), Pn a j Hn a ) x 



/l(Pl, Ml) •••/n (P )i[d Pl (2.7) 

where fj (jpj , fj,j ) are wave packets describing the momentum and spin distribution of the j th asymptotically bound 
cluster. 

The asymptotic Hilbert space, Jif, is the direct sum of all of the channel Hilbert spaces, including the one-cluster 
channels that correspond to iV-particle bound states, 

H f = ®iH ai , (2.8) 

and the multichannel injection operator $ : Hf — > H is 

$|f> IW (2-9) 

i 

where 

|f)=©o t |W- ( 2 - 10 ) 

For each partition a there may be 0, 1 or a finite number of channels. 

The unitary time evolution operator U ai (t) on each channel subspace, "H Qfc is 

U ak (t) = e-^P?/ 2 ^-^)* (2.11) 

where Pj is the total momentum of the j-th bound cluster in the channel , mj is the total mass of the j-th bound 
cluster of channel ctfc and — 6j is the binding energy of the j-th bound cluster of channel cti. 

The asymptotic time-evolution operator onJif is the direct sum of the channel time-evolution operators 

Uf{t)=®iU ai (t). (2.12) 
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The asymptotic Hamiltonian, Hf, is the infinitesimal generator of Uf(t). Multichannel M0ller wave operators 

n ± :H f ^H (2.13) 

are defined by the strong limits 

fi± = lim U(-t)<f>Uf(t) (2.14) 

t— f±CXD 

where U(t) = e~ lHt is the time evolution operator on H. The multichannel scattering operator S : Hf —> Hf, is 
defined by 

s = q\.Q— (2.15) 

In order to indicate the dependence of the wave operator f2± on H, Hf and $ I use the notation 

Q± = fl±(H,<f>,H f ) S(H,$,Hf) = Q,\_(H,<f>,Hf)Sl-(H,$,Hf). (2.16) 

I say that the scattering theory is asymptotically complete if the wave operators satisfy the following completeness 
relations: 

I n = n + (H, $, H f )Sl\.(H, H f ) = n_ (H, $, Hf){ll(H, $, H f ) (2.17) 

and 

I n , =nl(H,$,H f )n+(H,$,H f ) = nUH,$,H f )Sl-(H,$,H f ). (2.18) 

where In and In s are the identity operators on H and H / respectively. 
The intertwining relations, 

HQ±(H, $, H f ) = n ± (H, $, H f )Hf, (2.19) 
follow directly from the definition (|2.14[) and lead to energy conservation, 

[H f ,S]- = 0, (2.20) 

in the scattering operator. 

In all that follows the wave operators are assumed to exist and satisfy the completeness relations (|2.17[) and (|2.18l) . 
The two-Hilbert space formulation of multichannel scattering theory summarized above is equivalent to the standard 
formulation of multichannel scattering. It has the advantage that the notation allows all channels to be treated 
simultaneously. 

III. CLUSTER EXPANSIONS 

In this section I introduce combinatorial methods to treat cluster expansions in this work 000. These provide 
an efficient notation for computing the many-body interactions that appear in different equivalent Hamiltonians. 

I begin by introducing a useful notation. I let V denote the set of partitions of N particles into disjoint, non-empty 
clusters. I use lower case Latin letters, a, to denote partitions of N particles, n a to denote the number of clusters in 
the partition a, and n ai to denote the number of particles in the i th cluster of the partition a: 



N = J2n az . (3.1) 



!=1 

Thus a = (125)(37)(64) is a three-cluster partition of seven particles, with one three-particle cluster and two two- 
particle clusters. 

There is a natural partial ordering on the set of partitions of N particles given by 

a<Zb or b D a (3.2) 

if every particle in the same cluster of a is also in the same cluster of b. For example a = (125) (37) (64) C b = 
(125) (3467). 
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I let a U b denote the least upper bound of a and b with respect to this partial ordering and a H b denote the 
greatest lower bound of a and b with respect to this partial ordering. I let 1 denote the unique 1-cluster partition 
and denote the unique N cluster partition. For a = (125) (37) (64) and b — (125) (367) (4) these definitions imply 
a n b = (125) (37) (4) (6), a U b = (125)(3467), 1 = (1234567), and = (1)(2)(3)(4)(5)(6)(7). 

Next I introduce the operators that translate clusters. On each of the single-particle Hilbert-spaces, Hi, there is a 
trivial representation of the three-dimensional Euclidean group 

f/ i (x 1 i?)| P , At )= J2 \Rp^'y Rpx Di,,(R) (3.3) 

H'=—3 

where x are parameters of the space translation subgroup and R is a rotation. The matrix D 3 , (R) is the ordinary 
Wigncr function. I define 

£/ a (xi, Ri, ■ ■ ■ ,x„ a ,i?„ a ) := 

®n&i ®i 2 eo a U l2 (x 2 ,R 2 ) ■ ■ ■ ®j„ eo„ U ina (x na ,R na ). (3.4) 

These operators perform independent translations and rotations on the subsystems of particles in each cluster of the 
partition a. 

A bounded operator A on the ./V-particle Hilbert space has a cluster expansion if it can be expressed as a sum of 
terms associated with each partition a, 

A=Y^[A} a> (3.5) 

where the operators [A] a are invariant with respect to independent translations and rotations of the cluster of the 
partition a 

[L4] a ,E/ a (xi, ,x na ,Rn a )}- = 0, (3.6) 

and vanish when any pair of particles in the same cluster of the partition a are asymptotically separated: 

lim WWaUbixuRf- .x^.iZOMH =0 M a - (3-7) 

Equations (|3.6I and 13 . T[) provide a mathematical characterization of these two properties. When A is unbounded I 
will assume that these equations hold for a suitable dense set of vectors \ijj). 
For Da, Ub{- • • ) is a subgroup of U a {- • • ) so 

lim ||L4] a £/f,(xi, • • • ,x„ 6 ,i?„J|i/')|i = 
lim ||C^(x 1 „iJ 1 ... J x n6 ,iZ n J[A] ^)|| = ||[A]„^)|| b^a. (3.8) 

\Xi — Xj |— >oo 



It follows from (|3.6[) and (|3.7p that if A has a cluster expansion then 

lim IP-V^yt^X!,^,--- ,x„ b ,i?„J|^>|| =0 (3.9) 

Xi— X-j — »oo — 

which leads to the definition 

A b := (3.10) 

Do 

which is the part of A that is invariant with respect to translations of the individual clusters of 6, irrespective of the 
asymptotic properties: 

[A^Ub^Ru--- ,x nb ,R nb )}- = 0. (3.11) 
This is the part of A that remains after the clusters of the partition b are asymptotically separated. 
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It is also useful to define 

A b :=A-A b = Y}A]a (3.12) 

bjba 

which is the part of A that asymptotically vanishes when the different clusters of b are asymptotically separated: 

lim \\A b U b (x u R u --- ,x nb ,R nb )\i>)\\=Q. (3.13) 



|Xi— I— ^oo 



The incidence matrix, S a o b , called the zeta function of the partial ordering, a D b, has an inverse, called the Mobius 
function, S~^ b of the partial ordering, which also vanishes when a b. The Mobius function can be used to express 
[A] a in terms of A b : 



[A]a=J2 S aD b A b- ( 3 - 14 ) 



i _ / (-irnr=i (-r 6 - K-i)! «^ 



This inverse is explicitly known (Toj : 

^ = ( " -^"^ ^ (3-15) 

where n bi is the number of clusters of b contained in the i th cluster of a. 
The following identities are consequences of the definitions 

(AB) a = A a B a (AB) a = A a B a + A a B a + A a B a (A a ) b = A anb . (3.16) 

If A has a cluster expansion the connected part of A is the part of A that vanishes when any pair of particles is 
separated. It is 

[A] 1 = A-^2[A] a . (3.17) 

Using properties of the Mobius and zeta functions gives the following expression for [A]\: 

I4 = EZ s iDaSaD b [A] b = £ S^ a A a = S^A, + £ S^ a A a = (3.18) 

A + Y. 5 ^aAa (3.19) 

a#l 



where I have used the identities 



I define the coefficients 



from which it follows that 



= 1 A 1 =A. (3.20) 
C a :=-6^ a = (-r°(n a -l)\ (3.21) 



A=[A] 1 +Y J C a A a . (3.22) 



This separates the "connected" part [A]i of A from the disconnected part, Y^ajti^aA a , of A. 
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IV. SCATTERING EQUIVALENCES 

Of interest is a characterization of when two TV-body Hamiltonians are physically equivalent. It is customary in 
the literature to call two Hamiltonians equivalent if they are related by a unitary transformation. This is really 
insufficient. For a simple counter example consider two different short-ranged repulsive two-body interactions, V\ and 
V2 • The spectrum and multiplicity of the of the two-body Hamiltonians are identical. If the wave operators satisfy 
the completeness relations (|2.17H2.18"| then the operator 

A = Sl + (Hi,*,Ho)Cl\(H 2 ,*,Ho) (4.1) 

is a unitary operator on W. It follows from the intertwining relations, (|2.19[) , that A also satisfies AH2 = Hi A; 
however any two arbitrary repulsive potentials do not give the same phase shifts. So even though Hi and H2 are 
related by a unitary transformation, the scattering observables arc unrelated. For equivalent Hamiltonians I also need 
to require that the S matrix remains unchanged and the description of the free particles remains unchanged. 

The two-Hilbert formulation of scattering is useful in this regard. What is required in general is the unitary 
equivalence of the Hamiltonians 

H' = A^HA AA* = I (4.2) 

and S'-matrix equivalence 

S(H t $,Hf)=S(H' t & t H f ), (4.3) 

where Hf remains unchanged. Recall from the construction of the previous section that the operator $ also depends 
on H. 

To determine the requirements of S-matrix equivalence on A I use (| 2 - 161) in (|4.3|) to obtain 

tf + (H, H f )n_(H, $, H f ) = fll(H', H f )il- (H', H f ). (4.4) 
Using (|2.15|) in (|4.4|) gives the following candidate for A. 

A := Q + (H, $, H f )nl(H\ H f ) = (l_(H, $, H f )tf_{H', H f ). (4.5) 
The intertwining property (|2.19j) gives 

AH' = HA. (4.6) 
Unitarity of A follows from (|2.17H2.18| . which also can be used to show 

fl+(H,$,H f ) = {l + {H,*,H f )I nf = 

fi+(JT, $, H f )tf + (H', H f )il+(H', H f ) = 

Afl+(H',' <5>,H f ) = n + (AH'A\A$',H f ) =Q + (H,A$',Hf). (4.7) 
Subtracting the left from the right side of (|4.7j) and using the definition of the wave operators gives the identity 

= Sl+(H,*,Hf)-Sl + {H,A&,Hf) (4.8) 

which is equivalent to 

0= lim \\U(t)[$ - A&]U f (t)\f)\\ = Hm |p - A&]U f (t)\f)\\- (4.9) 

t— >oo t—>oo 

Similarly using the second equation (|4.4j) gives the corresponding relation with the other time limit 

0= lim || [$ - A<£>']Uf(t)\f}\\. (4.10) 

t— 00 

The vanishing of both time limits is important. The failure of S'-matrix equivalence in the case of the two repulsive 
potentials is because the two time limits lead to different unitary operators, A + 7^ A—, satisfying (|4.2[) . 



7 



The asymptotic conditions, (j4.9H4.10j) . along with the definition of the operator A (|4.5p . are consequences of the 
identity (|4.3p of the two scattering operators. 

Conversely, if both asymptotic conditions, (|4.9|) and (|4.10|) . hold for some unitary A then 

Q±(H,®,Hf) = Q±{H,A<f>',H f ) = AQ±(H',&,H f ). (4.11) 

Because this holds for the same A for both time limits it follows that 

S(H,$,H f ) = n\_(H,<f>,H f )Cl_(H,<£>,H f ) = 

n\(H' ,H f )A ] An_(H' \H f ) =S(H',&,H f ). (4.12) 

This shows that the asymptotic conditions (|4.9M.10]l are necessary and sufficient conditions for the invariance of the 
S- matrix. This result is the content of a theorem in formal scattering theory due to Ekstein (llj . 

I also need to determine the relation of $' to A and H in the context of Ekstein's theorem. I assume that A 
has a well-defined cluster expansion and I define an operator A a by turning off the parts of A that vanish when the 
clusters of the partition a are asymptotically separated. It follows that if I turn off the interactions between particles 
in different clusters of the partition a that H' will have the following limiting form 

H' = A^ HA -> H' a = A\H a A a (4.13) 

where H' a is a sum of transformed subsystem Hamiltonians associated with each cluster. It follows from the definitions 
(|2.7p and (|2.9[) that the channel injection operators $ Qi and associated with the partition a are related by 

= 4*<* (4-14) 



where $' Q . is an cigenstatc of H' a with eigenvalues 

E CH = 

Thus 



i=i 



$'=X)4 4 *o < . (4-16) 



Finally, if I want the subsystem Hamiltonians to be separately rotationally and translationally invariant, then each 
of the operators A a , obtained from A by turning off the parts of A that generate interactions between particles in 
different clusters of the partition a should also be translationally and rotationally invariant. 

Given a unitary A with a cluster expansion, equations (|4.14[) and (|4.16j) imply 

A^ = Y,AAi^ ai . (4.17) 

i 

The requirement that H and H' = A^ HA give the same S matrix is that A is a unitary transformation with a 
well-defined cluster expansion satisfying the asymptotic conditions 

0= lim \\Y][I ~ AAl]$ ai U ai (t)\f a M =0 (4.18) 

i 

for each channel, or equivalently because of the unitarity of A 

lim ||^[At-4]$ ai C/ ai (i)|f ai )||=0. (4.19) 

i 

If I assume that all of the f ? ; vanish except for the iV-body breakup channel then (|4=. 18[) implies 

0= lim || V[J- A^}U (t)\f )\\ =0. (4.20) 

i 

When A has a suitable cluster expansion, (|4.20p implies (|4.18p . This is discussed in section V. 
Equation (|4.20p is equivalent to 

0= lim || Y][I - A}U (t)\f }\\ =0. (4.21) 

i 

where we have used $o = I and Ho = H for the unique ^-cluster breakup channel. 

I refer to unitary transformations A satisfying (|4.21[) as scattering equivalences. It is easy to show that with this 
definition the set of scattering equivalences form a group with respect to operator multiplication. 
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V. ASYMPTOTIC PROPERTIES 



It is now possible to construct a parameterized set of scattering equivalences. Because the scattering equivalences 
A are unitary operators, it follows that A can be expressed as the Cayley transform of a Hermitian operator T 

a-££ r-n. (5.D 

In what follows I will assume that the Cayley transform, T, has a cluster expansion, 

r = £[r] a (5.2) 

where the [T] n are Hermitian, invariant with respect to translations and rotations of the clusters of a, and vanish 
when any of the particles in different clusters of a are asymptotically separated. Specifically 

lim miy/b^i?!--- ,x ni; ,i?„ b )|V)|| =0 b^a. (5.3) 

A sufficient condition to satisfy all of the cluster conditions is that [r] Q and T[T] a T, where T is the ./V-body kinetic 
energy operator, are both compact after one removes all of the momentum conserving delta functions. In what follows, 
rather than formally taking the cluster limit (|5.3[) . I use a switching parameter to turn off the parts of the operators 
that vanish in the cluster limit. Thus, to take the limit where the clusters of a partition a are separated, I formally 
write 

r(A) = r a + Ar a r(i) = r (5.4) 

and take the limit that A — > 0. I call this implementation of cluster properties algebraic clustering [l^; it separates 
the combinatorial aspects of cluster properties from the analytic aspects. 

While the cluster expansions are based on asymptotic properties of the operators A with respect to translations, 
the limits of interest in this paper are the time limits (|4.18H4.2lj) . Although I will not get into the technical details of 
the cluster limits, it is important to understand the relation between the cluster limit and the time limit. 

If I consider the time limit in equations (|4.19[) , it has the form 

lim - Al}<£ m U a M\U\\ = 

t-^-±oo * — * 



lim ||VA Q t e - 4E( <- ej) $ Qi |f Qi }|| =0 (5.5) 

t— »±oo *• — » 
i 

where the kinetic energy that appears in the exponent is the sum kinetic energies of each cluster of a. It looks similar 
to the cluster limit 

lim || VA a V Ep ^ Xj $ Q! |f Ql )j| =0. (5.6) 



To understand the relation between the limits in (|5.5j) and (|5.6|) I consider first a single degree of freedom. Consider 
the limit where the y-component of cluster i is being translated. The time limit above is bounded by a sum of terms 
of the form 



/+oo 



lim / f(p tv y^ x dp ty + lim / f(-p ty yP*y x dp ly . (5.7) 
A -*°°Jo x ->-°°Jo 

Both terms vanish by the Ricmann Lcbcsgue lemma if f(pi y ) and f(—Pi y ) are absolutely integrable on [0, oo]. 
The corresponding time limit (|4.21[) contains a term of the form 



A 



lim / f(pi y )e 2m ' X dp ly . (5.8) 
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If I let u = the time limit becomes 



lim / (/( V / 2^I) + /(-V2^))e 4tlA A /— du = lim / g(u)e iuX du (5.9) 

A->oo/__ V u X- 



wherc 



/ 777 ■ 

g{u):=(f(V2^) + f{-V*w))J— (5.10) 

V u 



is an absolutely intcgrable function of u if f(pi y ) is an absolutely integrable function pi y . Using an extension of this 
same argument it is possible to show that the two limits (|5.5|) and (|5.6[) arc equivalent, provided <£> ai , f ai and A a are 
all suitably well-behaved (i.e. so the resulting integrand is absolutely integrable). 

This means that the time limit associated with a given channel has a vanishing limit whenever the space limit 
associated with the same channel also vanishes. Once I eliminate the delta functions, the compactness condition 
always ensures that (|5.5[) and (|5.6p are satisfied. If the functions are smooth the fall-off is faster. 

For channels associated with the partition a the operators A a must vanish for both time limits. A sufficient condition 
for this to be satisfied for all partitions a is that 

lim \\{A-I)$ U {t)$ \f )\\=Q (5.11) 

t— ^±oo 

for the N cluster partition 0. In this case $o = I, Rq = R an d this condition becomes 

lim \\(A-I)U Q (t)$ Q \f )\\ =0 (5.12) 

t— f ±oo 

This ensures that I — A is a sum terms that vanish when all particles are asymptotically separated. This is the basis of 
our claim that (|4.21[) implies (|4.19[) and leads to the characterization (|4.2ip of the asymptotic properties of scattering 
equivalences. 



VI. CONSTRUCTION 



To construct a suitable class of operators A that can be used in variational calculations, consider operators A, 
where the Cayley transform has a cluster expansion 

A = \-§ r = rt - r = E^ (e- 1 ) 

aev 

where each [r] a is a Hcrmitian operator that commutes with U a {x.i, Ri, • • • ,x„ a ,i?„ a ) and satisfies the asymptotic 
condition 

lim H^^xi,,^--. ,x nb ,Rn b )M\\ =0 M«- (6-2) 

I also assume that after the momentum conserving delta functions are removed, the remainder is a compact operator 
with respect to the internal variables. This ensures that (|6.2[) holds. The means that the internal part has an expansion 
of the form 

[r]a = I x [f] a (6.3) 
where / is associated with the delta functions and compact remainder has the canonical form 

[r] a = ^|U)A n(U| (6.4) 

n 

with X an = A* n , limn^oo |A on | and (£ a m|£<m) = S mn . In a variational framework the coefficients X an and the 
orthogonal vectors \£ an } can be chosen to depend on variational parameters. 

In general the operators [T] , along with the original Hamiltonian, are the input to any calculation. In addition, 
because of Ekstcin's theorem, the transformation leads to scattering equivalent Hamiltonian characterized by a scat- 
tering equivalence with a cluster expansion of the form (|3.5[) . The Cayley transform may be unbounded, but it will 
have and algebraic cluster expansion of the above form. 
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The cluster expansion of T can be used to generate the cluster expansion of A. Since r a and [r] Q are related by 
the Mobius and zeta functions, it is possible to construct T a from the [r]&'s. 
I have 



r a = J2 5 <^bF}b (6.5) 



A a = (6.6) 
1 + iT a 

[A]a = J2 S ^b A »- ( 6J ) 

The AaS can be computed recursively on the number of clusters in the partition, starting with N— 1 cluster partitions. 

The nature of the general construction can be illustrated using a three-body example. In this case T(ij)(k) = [r](ij)(fc) ■ 
For the two cluster partitions, a = (ij)(k), I first solve the integral equation 

- — = - - i{ r ](ij)(k)- — • (6-8) 

For finite rank [T](y)(fc) this is an algebraic problem. [r] Q 's of the form (|6.4[) can be uniformly approximated by finite 
rank [T] a . For the special case that [r] a = AII a is a real constant multiplied by the direct product of the identity (in 
the conserved momentum variables) and a rank-one projection operator, equation (|6.8[) can be solved analytically. 
The solution is 

A (ij)(k) = TTTffl = TTTx 1X W)(*)- ( 6 - 9 > 

1 + l l L J(ij)(fc) i + lA 



To use these solutions to compute A I define 



R:= T+iT ^ 



R(ij)(k) — — (6.11) 



and 



For the special case of a rank one ^(ij)(k) 



/ + *[ r ](y)(fc) 

fl(l)(2)(3) =/■ (6-12) 



R m(k) - 1 - h , { -\ {k) — u m(k) ■ ( 6 - 13 ) 

1 + *^fe')(fc) 

In general the operators R^yA and R satisfy the resolvent identities 

R(ij)(k) = i - «[r](ij)(fc)-R(y)(fe) (6.14) 

and 

R = R(ij)(k) - iR (ij)(k)([r}(jk)(i) + [r](fei)(j) + [r](i23))-R (6.15) 

R = %)(2)(3) + (6.16) 
To get an equation for R note that (|3.15p and (|3 . 2 1 1) imply 

J2 Ca = L ( 6 - 17 ) 
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Using this with equations (13.21[) gives the following equation for R: 

R = j2c a R = J2 c « R "- i J2 CaR « TaR - ( 6 - 18 ) 

a^l a^l a^l 

Equation (|6 . is valid for any number of particles. For the three-particle case the driving term and kernel of (I6.18|) 
can be expressed in terms of the [r] a as 

^2C a R a = R(ij)(k) + R(ij)(k) + R(ij)(k) -21 = 

I - i[r](i2)(3)-R(12)(3) - *[r](23)(l)-R(23)(l) ~ «[r](31)(2)-R(31)(2) (6.19) 

and 

-iY,CaR a r a = 

a=£l 

i[r] (12 )(3)i?( 12 )(3))([r](23)(i) + [r](3i)(2) + [r](i 2 3)) 

- i[r](23)( 1 )i2(23)(l))([r](3i)( 2 ) + [r](i 2 )( 3 ) + [r](i23)) 
i[r]( 31 )(2)i?(3i)(2))([r]( 12 )(3) + [r] (23 )(i) + [r](i23)) 

+ *2/([r] (12)(3) + [r] (2 3) ( x)[r] (31)(2) + [r] (123) ) (6.20) 

= -*[r](l23) - [r](l2)(3)^(12)(3)([r](23)(l) + [ r ](31)(2) + [r](i23)) 
-[r]( 2 3)(l)-R(23)(l)([r](31)(2) + [ r ](12)(3) + [r](l23)) 

- [r](31)(2)-R(31)(2)([r](i2)(3) + I r ](23)(l) + [r](l23))- (6-21) 

The important observation is that this operator, which is the kernel of the integral equation (|6.18[) . is compact after 
delta functions that arise from overall translational invariance are removed. It follows that equation (|6.18[) can be 
solved my standard Frcdholm methods. The solution can then be used to construct A using 

A=\-^ = (l- l r)R= (6.22) 
1 + 11 

{I ~ i([T](l2)(3) + [r](23)(l) + I r ](31)(2) + [r](123))i?. (6.23) 

If the individual [r] a are finite rank (after all of the delta functions due to the translational symmetry are removed) 
then it follows that the kernel (|6.21[) is finite rank (after the overall momentum conserving delta function is removed) 
. This is most easily seen in the special case where all of the [r] a arc proportional to one-dimensional projectors (after 
the delta functions are removed). In this case the kernel (|6.2ip becomes 

-«A(123) 11(123) - (12)(3) II( 1 2)(3))(A(23)(i)II(23)(i) + A( 31 )(2)II(3 1) ( 2 ) A( 123 ) il( 12 3) 

1 + ZA(i2)(3) 

n(23)(l))(A(31)( 2 )II( 31 )( 2 ) + A( 1 2)(3)II( 12 )(3)A( 1 23)II( 1 23) 



1 + *A(23)(i) 
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-n(31)(2))(A(12)(3)II(i2)(3) + A(23) (1) 11(23) (1)A(123) 11(123). ( 6 -24) 



1 + *A(3i)(2) 

This is a finite dimensional matrix involving the ten operators n (123 ), n( 12 )(3)II(23)(i), n^)^)!!^!)^), n( 12 )(3)n( 12 3), 
n(23)(i)II(3i)(2), n(23)(i)II( 12 )(3), n( 23 )(i)n( 12 3) , n ( 3 1 )( 2 )II( 2 3)( 1 ), n (31 )( 2 )n( 12 )( 3 ), and n( 31 )( 2 )IT( 12 3). After the overall 
momentum conserving delta function is removed the range of the this operator is a ten dimensional vector space. The 
resulting integral equation involves solving a system of 10 linear equations. When the operators [r] a are finite rank, 
rather than rank one, the matrix is larger, but it is still finite dimensional. 

This construction can be extended to any number of particles. The kernel of the integral equation for the iV-body 
R is still finite rank if all of the input [r] a are finite rank. Thus, for finite rank [T] a the construction of A involves 
only quadratures and linear algebra. 

Returning to the three-body example have 

H = T + V(i2)(3) + V(23)(l) + V(3i)(2) + V(i 23 ) (6.25) 

Ai)(2)(3) = / (6-26) 

[ J 4](12)(3) = (I - i[T](l2)(3))-R(12)(3) - 1 = - 2 *'[r](i 2 )(3)-R(12)(3) (6.27) 
[^](23)(1) = (I - ?; [r](23)(l))i?(23)(l) - I = -2i[r]( 23 )(i)J?(23)(l) (6.28) 
[^](31)(2) = (I - «[r]( 3 l)(2))i?( 3 l)(2) - I = -2i[r]( 31 )( 2) i?(3i ) (2) (6.29) 

[4(123) = -lJ2 C « R " TaR = ( 6 ' 3 °) 

#(12) (3) (PI (23) (1) + I r ](31)(2) + [ r ](123))-R+ 
-R(23)(l)([r](31)(2) + [r](12)(3) + [r](l23))-R+ 
^(31)(2)([r](12)( 3 ) + [r](23)(l) + [r](123))-R 

- 2([r] (12) (3) + [r] ( 23)(i) + [r] (31 ) ( 2 ) + [r] (123) )i? (6.31) 

I can use these cluster expansions of the scattering equivalences to determine the cluster expansion of the transformed 
three-body Hamiltonian 

H' = A^ HA (6.32) 
H' a = A\H a A a (6.33) 

i H 'h =H' - Y, C - H 'a = A^ HA -Y^C a AiH a A a . (6.34) 

a/l a^l 

This means that the transformed two-body interactions are 

V( ijm = [H% m) = (I+ [A]\ mj) )(T + %) W )(I + [A] (ijm ) - T = (6.35) 
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l A ]\ij)u)) v (,mk) i A hjm + %)(fe) [ A hjm + T \ A hj)U) ( 6 - 36 ) 

where T is the three-body kinetic energy and the A^^) arc given by (|6.27ti6.29|) . 

An important observation is that VL.^,^ only depends on T, V^^, and [T]^j^ k y It docs not depend on [r]( 123 ). 
This means that after one chooses [rjfo'Vfc) to give a transformed two-body interaction, it is still possible to use 
the freedom to independently choose [r]( 123 ) to transform the resulting three-body interaction without changing the 
transformed two-body interactions. 

The transformed three-body interaction is 

V(W = A ^ HA " F (23)(l) - ^(31)(2) = [AHA]! = (6.37) 

[A] {123) HA^ + 

[4(12)(3) V(12)(3))([^](23)(1) + [4(31)(2) + [^(123)) + (V(a3)(l) + V(31)(2) + %23))^ 

[4(23)(1) ((T + V (2 3)(1))(L4]|3 1)(2) + [A][ 12)(3) + [A]| 123) ) + (F ( 31) (2 ) + V (12)(3 ) + ^(123))^ 

[4(3i)(2) ((T + F ( 3i)(2))([A]J 12)(3) + [A]| 23)(1) + [A]J 123) ) + (F (12)(3) + V {23)(1) + V (m )A^ + 
T[A] 

(123) + %2)(3)([4(23)(1) + [4(31)(2) + [4(123)) + 
^(23)(l)([4j 3 l)(2 ) + [4<12)(3) + [4(123)) + 

^(31)(2)([4fl2 )(3) + WUW) + [4(123))- ( 6 - 38 ) 

This is expressed as a sum of completely connected terms; it could be expressed in a more symmetric form but 
that would involve more terms. The entire expression depends on the operators [r] a that depend on the variational 
parameters. 

If the [-^](ij)(fe) have already been determined by fixing the two-body interaction then one can start with the 
transformed potential and use an A where only [r]( 123 )] is non-zero to get an optimized three-body interaction. 
Alternatively one could state with the original potential and leave the [A]^)^) fixed in the above expression, with 
all of the variational parameters in [r]( 123 ). 

VII. CONTROLLING THE HAMILTONIAN 

I order to use variational methods to determine the best choice of Hamiltonian a positive functional is needed that 
can be minimized. It is possible to either work recursively on the number of particles, by first determining two-body 
interactions, followed by the three-body interaction, or alternatively to determine all interactions simultaneously. 

The simplest type of functionals arc of the form 

F(V) = Tr(pV^V) 1 / 2 (7.1) 

where p is a positive, rotationally and translationally invariant operator. The trace is only taken over the variables 
that remain after the momentum conserving delta functions are removed. Thus for two-body interactions of the form 

V(k,k', m --- m )6(p' -p) (7.2) 

F(V) would have the general form 

F(V) = J dkdk'dk'V(k,k',7 /1 .-. w )F*(k",k',7 ?1 ...7 ?A r)p(k,k") (7.3) 
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with obvious generalizations for three-body interaction. 



F(V) = J dkdk'dk"V(k - k', m ■ ■ ■ m)V*{k! - k", • • • r] N )p(k" - k') (7.4) 

If the starting potential is local this expression has to be modified because VV^ is a function of the difference k' k 
which leads to an infinite volume factor. While local potentials can be treated by using a different positive functional, 
an alternative is to note that if V = Vi oc + V r then F(V) = F{Vi oc ) + F'(V r , Vi oc ). It is only the first term that is 
infinite, but this term does not depend on the variational parameters. The second term will be finite for suitable A 
and it contains all of the dependence on the variational parameters. It follows that the critical value of the variational 
parameters can be determined by requiring that all partial derivatives of the second term at the critical value of the 
parameters. 

The general procedure is to start from a given iV-body Hamiltonian, H , and a parameterized set of scattering 
equivalences A(r/i, • • • , r) n ) where r\i are variational parameters. The scattering equivalences A(r)x, • • • , rj n ) generate a 
parameterized set of equivalent Hamiltonians: 

H 'ivi, ■ ■ ■ ,Vn) = -4 T (?7i, • • • ,r) n )HA(rn, ■ ■ ■ ,r?„). (7.5) 

They have cluster expansions 

H> = T + Y,V! j {r h ,--- ,Vn) + J2 V m(riu--- ,Vn) + ■ • • + Vjffr, ■ ■ ■ , r,n)- (7.6) 

ij ijk 

The two, three, four • • • N-body interactions all depend on the choice of variational parameters. 

For example, to construct two-body interactions that have primarily low-momentum content I would choose a 
functional that is large when the momenta are large. The functional has to be chosen so the trace is finite for all 
interactions in the model space. 

A functional of the form 

p(k, k') = tanh(a + k 2 /ko)tanh(a + k' 2 /ko), (7.7) 

where a is a small dimensionless quantity, would suppress momentum components above the scale kg. Alternatively 
I can design positive functional that weaken three-body forces or reduce two-body correlations. 
Finding minimum of the functional 

Tr(pVl 2 (m, ■ ■ ■ ,Vn)V 12 (m, ■ ■ ■ ,r)n)) (7.8) 

with respect to the parameters 771, •• • , rj n selects equivalent potentials that have low-momentum content. 

After the two, three, ■ • • N— body interactions have been determined, then I can use the new Hamiltonian as the 
starting point. I can construct a new set of interactions using scattering equivalences with [r]^^) = 0. These 
scattering equivalences only affect the three and more-body interactions. I can choose a new three-body p that 
emphasizes some desirable feature of the three-body interaction. The local minimum generates a new three-body 
interaction. Combining the two scattering equivalences leads to an scattering equivalence A that transforms H" = 
A^HA. 

If this is embedded in the A^-particle Hilbert space it (1) generates the selected two and three-body interactions, (2) 
new 4, 5 • • • iV-body interactions, and (3) explicit unitary transformations, A, that can be used to generate transformed 
operators like electromagnetic current operators 

J»'(x) = AW x {x)A. (7.9) 

VIII. SIMPLE EXAMPLE 

To illustrate the method I consider a two-body Hamiltonian of the form 

H=^ + V t (8.1) 
where I assume that V is a local potential. I consider a parameterized rank one unitary transformation of the form 

A(X) = /+|.g) 1 _^ g|g) (.9| =I + \g)f(\)(g\ (8.2) 
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where \g) is a fixed form factor and A is a variational parameter. The transformed potential is 

V'{\) = A\\)HA{\)-¥-. (8.3) 

The transformed potential differs from the original potential by the addition of a finite number of separable terms. It 
has the form 

V'(X) = V + \g)f*(\)(g\H + H\g)f(X)(g\ + \g)r(X)(g\H\g)f(X)(g\ =V + V r (X) (8.4) 

The first term in this expression is local but independent of A. The remaining terms are separable and depend on A. 
I use the density (|7.7p , with a chosen value of ko . It has the form 

P=\X)(X\ (8.5) 

leads to the variational function 

F(X) := ( X \ (Ft'(A)F'(A) - V*V) \x) = (x| (W r (A) + V}(X)V + (X)V r (X)) \ X ). (8.6) 

The subtracted contribution, V^V, eliminates the infinite constant that appears for local V. The terms in the resulting 
expression are 

F(X) = 

(x\V^\g).r(X)(g\H\ X ) + (x\V^H\g)f(X)(g\ X ) + ( X \V^\g)r(X)(g\H\g)f(X)(g\ X ) + 
( X \g)r(\)(g\HV\x) + (x\H\g)f(X)(g\V\ X ) + (x\g)f* (X)(g\H\g)f(X)(g\V\x) + 
({ X \g)f*(X)(g\H + ( X \H\g)f(X)(g\ + ( x \g) f (X){g\H\g) f(X)(g\) x 

(\g)r(X)(g\H\ x ) + H\g)f(X){g\ X ) + \g)r(X)(g\H\g)f(X)(g\ X )) . (8.7) 

This has the form 

F(X) = d/(A) + c*J*(X) + ca/(A)/*(A) + c 3 / 2 (A)/*(A) + C ;/* 2 (A)/(A) + c 4 (/(A)/*(A)) 2 (8.8) 

with 

The coefficients Ck are liner combinations of the integrals (xl^lff)i (g\H\x), (x\^H\g), (g\HV\x), (g\x), (g\H\g), 
(g\g) and (g\H 2 \g). Since these do not involve A they only have to be computed once. Although /(A) is complex, 
F(X) is a real function of A. The A dependence is a rational function. 

The critical value of A = A c is determined by solving ^j(A c ) = for A c . The resulting transformed Hamiltonian 

H'=^-+V'(X C ) (8.10) 

gives the same binding energies and phase shifts as the original potential of any value of A. The critical value of A 
will lead to a potential that suppress momenta above kp. Obviously a softer potential will result if a larger class of 
unitary transformations A are used. 

The original Hamiltonian did not have to be diagonalizcd to find the new potential. In this case, by varying A from 
to its critical value it is possible to continuously evolve the initial local potential to the final soft potential. 

Since the unitary scattering equivalence is given as an explicit operator valued function of A, I can calculate how 
observables evolve with the parameter A. For example the electromagnetic current operators transforms as follows: 



J»(x)' = J^(x) + r(X)\g)(g\J^x) + f(X)J»(x)\g) (g\ + f* (X)f(X)\g)(g\J^x)\g)(g\ (8.11) 
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Finally, give the two-body unitary transformation for each pair of particles, Ay (A), it is possible to construct the 
corresponding three-body unitary operator following the method of the previous section. In terms of the above 
parameters, for three identical particles A has the form 



with 



2 + f(X)(g\g) 



I — ioL , 
A = —— 8.12 
I + ict 



(|012)(ffl2| + |523>(323| + |S3l)(ff3l|) ■ (8.13) 



If the symmetric product of this unitary transformation for each pair is applied to the corresponding three-body 
Hamiltonian the transformed three-body Hamiltonian will have the form 

H' = A\\)HA = K + y/ 2 (A) + 143(A) + V^X) + V{ 23 (\). (8.14) 

The three-body force terms will appear even if the original Hamiltonian has only two-body forces. The computation 
of A from (|8.12[) involves quadratures and linear algebra, as discussed in section VI. 

The evolution of the current and the three-body Hamiltonian from their original to their final values can be 
determined by varying A from zero to the critical value, Aq. 



IX. CONCLUSION 



In this paper I determined conditions that are necessary and sufficient for two Hamiltonains to be physically 
equivalent. I used the characterization of these unitary operators to construct a large class of equivalent TV-body 
Hamiltonians that depend on variational parameters. There is considerable freedom in choosing the space of equivalent 
Hamiltonians. By choosing functions whose local minima select Hamiltonians with desirable properties from the 
space of equivalent Hamiltonians, it is possible to select classes of equivalent potentials with desirable properties. The 
general freedom available allows for the possibility of selecting two-body interactions with desirable properties, then 
subsequently selecting among equivalent three-body interactions with desirable properties. This procedure can be 
continued for any number of particles, allowing independent control of the two, three, four, • • • interactions. Because 
the fc-body parts of A affect all operators with k or more particles, one hopes that desirable properties of the fc-body 
interaction might persist for the k + m body problems. 

While in general it is possible to systematically weaken three and more-body interactions using these methods, it 
is not generally possible to eliminate them. The extremal interactions that are generated are not fundamental, they 
depend specifically on the choice of positive functional that is used to select these interactions. 

The selection of equivalent potentials does not require diagonalizing any Hamiltonians; it only requires finding 
local minima of some user defined functionals. The functionals are designed so they get large for interactions with 
undesirable features. Once the operators [r] a are determined variationally, it is then possible to construct scattering 
equivalences A that operate on systems of any number of particles, and can be used to construct equivalent observables 
in the transformed representation. For a large class of variational Hamiltonians the operators A can be constructed 
from the r a by finite linear algebra. 

The general method can be combined with other methods, such are renormalization group methods, to reduce the 
strength of the transformed three-body force without changing the transformed two-body interactions. 

The characterization of the group of scattering equivalences demonstrates the large class of equivalent Hamiltonains 
that can be selected by considering only spectral properties and scattering observables. This leads to a lot of flexibility 
in building equivalent models of the quantum A^-body problem. 

The author would like to acknowledge support for this work by the U. S. Department of Energy, Office of Nuclear 
Physics, under contract No. DE-FG02-86ER40286. 
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